Anisotropy due to Brans-Dicke Theory 
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It is well known that the viscous Bianchi type-I metric of the Kasner form is not able to describe 
an anisotropic universe, which satisfies the second law of thermodynamics and the dominant energy 
condition in Einstein's theory of gravity. We investigate this problem in Brans-Dicke theory of 
gravity using the Bianchi type-I metric with the perfect fluid. We show that it is possible to have 
the dominant energy condition and the growth of entropy in this model. Also we apply this model 
to explain the anomaly concerning the low quadrupole amplitude of the angular power spectrum of 
the temperature anisotropy observed in WMAP data. 
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One of the most interesting and complete scalar tensor theory of gravities is Brans-Dicke (BD) theory depending 
on the BD parameter wbd and the gravitational constant G is replaced with the reciprocal of a scalar field <I> [l| . 

Based on the Einstein's gravity theory, it is known that a viscous Bianchi type-I metric of the Kasner form is not 
possible to describe an anisotropic universe model which satisfies both the dominant energy condition (DEC) and the 
second law of thermodynamics [3, [1] . This problem might be resolved in some specific scalar-tensor theories of gravity 

The one of possible applications of the anisotropic universe to the phenomena is the explanation of anomalous 
feature of low quadrupole moment of the cosmic microwave background (CMB) anisotropy. If the large scale spatial 
geometry of the universe has the plane symmetric with an eccentricity at the last scattering of order 10~^, then the 
quadrupole amplitude of the CMB can be reduced with respect to the value of the best fit ACDM model without 
affecting higher multipoles of the angular power spectrum of the temperature anisotropy 0, |M 0] • 

The Bianchi type-I solutions of BD theory including perfect fluid is also well studied [3, Si- In what follows, we 
analyze in detail the solutions of this model and investigate both the entropy growth and the DEC of it. 

This paper is organized as follows. In the next section we investigate the solution of BD theory in a Bianchi type-I 
metric with perfect fluid. We analyze the various energy conditions and the growth of entropy of this model in section 
3. We also consider the phenomena and the possible application of this model to the low quadrupole problem of the 
CMB in section 4. In section 5 we reach our conclusions. 
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2. BRANS-DICKE FIELD EQUATIONS 

We start with the Bianchi type-I metric where the line element is given by 

ds^ = -c^dt'^ + a,{tfdxl 



(2.1) 



i=l,2,3 

where ai is the scale factor of each spatial direction. BD theory is described by the action 

5 = /" d^a 



(2.2) 



where $ is a scalar field, wbd is the BD parameter, and >Cfluid is the Lagrangian of the ordinary matter component. 
From the above action (|2.2[) . we can find the BD field equations 
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The field equations (|2.3p are explicitly given in terms of the metric 
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where ^ = 010203. From the Bianchi identity we get 

P'^V = p-i^^ni , (2.9) 

where pj^j and T^ni are the matter density and the volume element at a given time ^inj. 
The shear tensor has the form 

<y^.u = hlU^^.^pX-\^^^^-^ (2-10) 

where is the four velocity, which satisfies U^Ui_i = —1, 6 = Uf"^ is the scalar expansion where the semicolon means 
the covariant derivative, and h^^ is the projection tensor defined as h^^i, = g^^i, + Uf^U^. We can identify the scalar 
expansion 9 as the mean Hubble expansion rate H by 

3 \ fli 02 03 y iv 

Now we can find the shear tensors from the above definition (j2.10p 

au = \M^ + ^-2^\ cToo-0. (2.12) 

Later we will investigate the effect of anisotropy expansion in the CMB temperature anisotropy. Thus we are only 
interested in the matter dominated epoch in this model. In this case, we can choose the initial time iini as icq when 
the matter and the radiation components have the equal amount of the energy density. Thus, we can find the value 
of tini from the current observation from the Wilkinson Microwave Anisotropy Probe (WMAP) Q 
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It is well known that the above field equations have analytic solution in this epoch ff, . Equations (|2.8|) and (|2.9 
give the relations 
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By integrating equation (|2.7|) with taking into account the Eqs. (|2.14|) and (|2.15p . we find 

l/(3cJBD+4) 



r(i2 + 2Tt) 



(2.14) 
(2.15) 

(2.16) 



where r = (47rp?^;/$i,ii)[(3a;BD + 4:)/{2ujbu + 3)] and T is an arbitrary constant with the dimension of time. We 
can find the limit on the value of T in section 4 in order to be consistent with data. The metric coefficients Ui are 
evaluated from the equations (|2J| . (|2T4l) . ((2T5l) . and ((2T6l) to give 



a,{t) = ai{t ini) 
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(2.17) 
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where ai(tini) and are arbitrary constants satisfying 

ai(tini)a2(iini)a3(iini) = V^ni = 1 , Cl + 6 + ?3 



0. 



(2.18) 



From the above equation (|2.17p . we can understand the isotropy of the universe. When the age of the universe is 
about the same order as T, then we should have the small ^ in order not to have the huge anisotropy of the expansion. 
However, as the universe gets older, this contribution becomes smaller and we will see the almost isotropic universe 
independent of the value of ^. 

We can derive the energy density and the pressure of BD field from the Eq. (|2.5p 
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From these equations (|2.19p and (|2.20p , we can define the directional components of the equation of state (eos) of BD 
field 
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It will be useful to define the sum of the pressure parameters by 
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FIG. 1: The directional components (Ais) and sum of them (A) of the eos of BD for the different values of T when we use 
choose the value of ubd as 120. a) The directional components of the eos of BD field (Ai) for the two specific values of T. 
When T = 3.12 x 10"^°to, all of AiS converge to the one value —0.4 during the interested epoch, b) The sum of Ai indeed 
converge to -1.2 when T = 3.12 x 10"^°to. 



In Fig[TJ we show the evolutions of A^s and A for wbd = 120 which is chosen by the CMB power spectrum constrain 
for ujbd (> 120) in the matter dominated epoch [lo[. We choose to as the present age of the universe and T as 
3.1210^^"io or 2.5 x 10~^to in this figure. The first choice of T value is given to be consistent with data and the 
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second one is for the demonstration for T with the similar magnitude as tinj. We also specify the planar symmetry 
universe (i.e. ai = a2 = a and as = b) in this figure, which constrains the value of ^^s as shown in the next section. 
The directional components of the eos of BD field (A^s) converge to one finite value (—0.4) independent with the value 
of wbd when the time increases. Thus, the sum of each component of eos which is identical to the eos of BD field 
(A) also converges to the finite value (—1.2) near the present. These are shown in the above figure. In the left panel 
of the figure [1] the solid and dotted lines show the evolutions of Ai (same as A2) and A3 when T = 2.55 x 10~^to. 
In this case we can see the anisotropic pressure contributions on each direction. However, dashed line indicates the 
both Ai and A3 when T ~ 3.12 x lO^^^io- They show the identical evolution in this case and the anisotropy of the 
universe is very small. The right panel of the above figure shows the evolution of A for the different values of T. 
We can find two interesting features from this figure. The first, in the matter dominated epoch all of the pressure 
contribution comes only from the BD field because the pressure of matter is zero. So anisotropy will be decreased 
because each directional dependent pressure components of BD field converge to the same value. The other is that 
the energy density of BD field is increased as the universe expanded {p^^ oc a°-^) instead of decreasing because the 
effective eos of BD field is smaller than —1. 



ENERGY CONDITIONS AND THERMODYNAMICS 



The energy-momentum conservation of total energy (V^ T^^/^+T/f = 0) is a consistency condition originating 



from the geometric Bianchi identity (V^G(^ = 0) with being the Einstein tensor given in Eq. 
equations ([^3]) . (|2.19p . and (|2.20p . then we have 
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Now we can express several useful quantities explicitly by using the equations from p.l4p to p.l7p 
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After we use the above equations (|3.2I 
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(|3.5p into the Bianchi identity p.ip . we have 
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To satisfy the above equation, we have 
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Again when we adopt the CMB power spectrum constrain for o^bd (> 120) in the matter dominated epoch, then we 
can abort the first condition in the above equation (|3.7p . If we consider the ellipsoidal universe {i.e. ai = 02 = a and 
03 = 6), then we can find the values of S^i from the equations (|2.18p and 
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FIG. 2: The values of coefficients of equations (|3.13|l and (|3.15p when we use the constraint ljbd > 120. a) To get the positive 
energy density, the investigated time should be bigger than 0.0034T. b) The SEC can be satisfied when the observing time is 
bigger than 0.0055T. 



The total energy density, the pressures, and the sum of pressures can be explicitly represented by 
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From the above equations (j3.10p - (|3.12p we can find the constraints of t to get the positive definite of energy density, 
dominant energy condition (DEC), and strong energy condition (SEC) 
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As we can see in the left panel of the figure [H if we restrict wbd > 120, then we can satisfy the positive energy 
density constraint when the observing time is greater than 0.0034T. The DEC is always satisfied for any time interval. 
The SEC is contented with t > 0.0055r. 

In general the energy momentum tensor for a viscous fluid is given by 



T — 



p+ {p- eO) UfJJy + (p - ee)gf,y ~ 2'qaf,y 



(3.16) 



where e and rj arc the bulk and shear viscosities, respectively. If we use the equation (|2.20p and define Pf^ — pi — 2rja\ 
(no sum over i) then we can find the shear viscosity as 
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The second law of thermodynamics requires 77 > and this condition is satisfied by given ^s p.9p and T = 3.12 x 
lO^^^io- Thus this model satisfies both various energy conditions and the second law of thermodynamics. 
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4. PHENOMENA 




FIG. 3: The time evolution of the ratio ai{t) to aj{t) for the different values of i^bd with T = 3.12 x lO^^^to. The anisotropies 
on each direction get smaller as the universe evolves, a) The evolution of ai{t) /ai{t) for ojbd ~ 120 (solid line) and 200 (dashed 
line), b) a3(t)/ai(t) for the same values of cjbd with same notation as in the left panel. 



In this section we investigate the cosmological models with planar symmetry. In this case, we have the relation 
between ^^s as given by the equation p.9p . If we check the scale factor in this epoch, then from the equation (|2.17p 
we have 
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where r = t/^o, To = T f^^id is again the present age of the universe. If the universe is isotropic and matter 
dominated, then the scale factor evolves as a/(i) = a/(io)(V^o)'^^^- H we set a/(io) = ai(to) = Ij then we can show 
the ratio between (t) and a/ (t) as in the figure [31 We show the ratio of ai {t) jaj {£) in the left panel of the figure 
[3l We choose Tq = 3.12 x 10"^" in this figure same as in the previous section. As we expected, the anisotropy of the 
universe is getting smaller as time evolves and it will become isotropic at near present. For the bigger values of wbd, 
we have smaller anisotropies. In the right panel of the figure [3l we also show the ratio of 03 {£) / ai (t) . 
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FIG. 4: The cosmological evolution of ratio to p™/"!? for the different values of lobds with T — 2.55 x 10 ^to and 
3.12 X 10~^°io- When T = 3.12 x 10"^°to, the value of /5^°<l>/p™ converges to 5/(6i.jbd) as we can find in the equation (|42)) . 
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If we check the evolution of the energy density of matter and BD field, then from the equations (|2.15p and (|2.19p 
we have the absolute values of the ratio of them 
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As we can see in the figure [H the absolute value of the energy density of BD field is decreased as time goes compared 
to that of matter. From the equation (|2.17p . we can see that as time goes the universe will be more isotropic and we 
will hardly see anisotropy. If the matter field keeps the dominant component and the universe becomes isotropic, then 
we have the matter dominated epoch with scaling a ~ t^^^. The Brans-Dicke field dominated epoch in the isotropic 
universe is investigated intensively in the literature In this case, the universe undergoes zero acceleration epoch 
a ^ t. 

If we assume that the universe is plane symmetric with an eccentricity at the last scattering surface (edoc) of 
order 10^^, then the quadrupole amplitude of the CMB temperature fluctuation can be reduced with respect to 
the observed value of the WMAP without affecting the amplitudes of higher multipoles fi|. The temperature 
anisotropy measured in a given direction of the sky can be expanded in spherical harmonics as 

= Va/mi1m("), (4.3) 

where fi is the direction of the photon momentum, Tg ~ 2.TiK is the present average temperature of CMB radiation, 
and aim are the multipoles. From this, the power spectrum is given by 

(4.4) 

In the plane symmetric universe with a small eccentricity, the temperature anisotropy is a linear superposition of the 
isotropic temperature fluctuation AT/ and fluctuations due to the anisotropic background AT^ 

AT = AT/ + ATa . (4.5) 

We may write the amplitude 

aim = aL + ■ (4-6) 

From the null geodesic equation, a photon emitted at the last scattering surface with temperature Tdcc reaches to the 
observer with the temperature 

To(n)=ro(l-ieL^2)^ (47) 




where Tq = Tdocfldoc and edoc — y (ff^) ~ ^- Therefore, the temperature anisotropy due to the anisotropic back- 
ground is given by 

AT^ 1,9 , ^ 

-jA = 2^doc«3 ■ (4.8) 

We need to get edcc — 0.67 x 10~^ in order to get the proper mean value of observed quadrupole anisotropy 



In the figure [5] we show the cosmological evolution of the eccentricity of the model (e) for the ujbu = 120. As we 
show in the figure when T = 3.12 x 10^^°to we can obtain the expected value of edoc — 0.67 x 10^^. 
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FIG. 5: The cosmological evolution of e{t) for T = 3.12 x 10~^°to with a;BD = 120. 



5. CONCLUSION 



We show that the Bianchi type-I metric in Brans-Dicke theory of gravity with matter fluid can satisfy the second 
law of thermodynamics and the various energy conditions. The cosmological evolution of the Brans-Dicke field reduces 
the possible primordial anisotropy of the universe. 

We have considered the alternative explanation for the anomaly in the quadrupole amplitude of CMB temperature 
spectrum. A small eccentricity due to the Brans-Dicke field may solve this conundrum. 
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